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ABSTRACT 

We obtained four Chandra/ ACIS-S observations beginning two weeks after 
the end of the November 2000 outburst of the neutron star (NS) transient 
Aql X— 1. Over the five month span in quiescence, the X-ray spectra are 
consistent with thermal emission from a NS with a pure hydrogen photosphere 
and R 00 =lh.Q 4 L%% (d/5 kpc) km at the optically implied X-ray column density. 
We also detect a hard power-law tail during two of the four observations. The 
intensity of Aql X— 1 first decreased by 50±4% over three months, then increased 
by 35±5% in one month, and then remained constant (<6% change) over the 
last month. These variations in the first two observations cannot be explained 
by a change in the power-law spectral component, nor in the X-ray column 
density. Presuming that Roo is not variable and a pure hydrogen atmosphere, 
the long-term changes can only be explained by variations in the NS effective 
temperature, from /cT e fj )00 = 1301 3 ; eV, down to 113±4 eV, finally increasing to 
118^4 eV for the final two observations. During one of these observations, 
we observe two phenomena which were previously suggested as indicators of 
quiescent accretion onto the NS: short-timescale (<10 4 sec) variability (at 32+g% 
rms), and a possible absorption feature near 0.5 keV. The possible absorption 
feature can potentially be explained as due to a time-variable response in 
the ACIS detector. Even so, such a feature has not been detected previously 
from a NS, and if confirmed and identified, can be exploited for simultaneous 
measurements of the photospheric redshift and NS radius. 
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1. Introduction 



Brown, Bildsten & Rutledge ( |1998| , BBR98 hereafter) showed that the core of a 
transiently accreting neutron star (NS), such as Aql X-l (for reviews of transient neutron 
stars, see |(Jhen et al. 1997] ; |Uampana et al. 1998a| ), is heated to a steady-state temperature 



by nuclear reactions deep in the crust during the accretion outbursts. The timescale for the 
core to reach this steady-state is ~10 4 yr (see also |Cofpi et al. 2000| ), after which the NS 
emits a thermal luminosity in quiescence of (BBR98) 



(M) \ Q_ 

10- l0 M Q jY- l J 1.45MeV/m f 



L q = 8.7 x 10 33 ( 1M _ in X 7/._„_ 1 ) T^TTTT— ergs s" 



where (M) is the time-averaged mass- accretion rate onto the NS, and Q is the amount of 
heat deposited in the crust per accreted nucleon ( [Haensel fc Zdunik 199C| see [Bildsten 



Rutledge 2000| for a discussion). Rutledge et al. (|2001a| , R01 hereafter) showed that, in the 



case of Aql X— 1, the quiescent X-ray luminosity was within observational uncertainties of 
that predicted on the basis of its time-averaged accretion luminosity. This observation, and 
similar observations of other quiescent transient neutron stars (qNSs), strongly supports 
the scenario that deep crustal heating provides a "rock bottom" quiescent luminosity. 

This quiescent thermal emission does not preclude that accretion onto the NS surface 
(at a rate M c ) continues in quiescence at a substantially lower rate than during outburst. 
[Narayan et al. (1997)| argued that the analogous black hole (BH) systems accrete during 



quiescence via an advective flow that radiates inefficiently and makes quiescent BHs dim 
or non-detectable. This same accretion flow onto the NSs should be observable, with 
efficiencies e = 0.20 (Lx ~ eM c c 2 ). Statistical comparisons show that the BHs are, on 
average, less luminous than NSs by a factor of ~10 ( |B arret fc Vedrenne 1994| ; |IN arayan] 



|et al. 19971 ; |Asai et al. 1998|) , although in individual cases they may be more luminous (V404 
Cyg is an example of a BH which is more luminous than most qNSs; [Kong et al. 2002|) . 
While this average difference in luminosity argues for a difference between the emission 
properties of NSs and of BHs, perhaps due to the presence of a surface in the former, 
detailed X-ray spectroscopy is still required to distinguish quiescent NSs from quiescent 
BHs in individual cases where the compact object mass is not measured nor type-I X-ray 
bursts observed ( Rutledge et al. 2000| ). 
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Assessing the relevance of accretion as a source of the NS quiescent emission is 
hampered by the difficulty of predicting M c . In the context of an advection-dominated 
accretion flow (ADAF; [INarayan fc Yi 19"9"3| ), M c is set by the evaporation rate from the disk 
which is in turn related to the unknown disk truncation radius. Similar limitations apply for 
other low M solutions, such as the Adiabatic Inflow/Outflow Solution (ADIOS; [Blandford 



fc Begelman 1999|) . We are thus searching for direct indicators of quiescent accretion onto 



the NS: short term (timescales much less than the orbital period) intensity variability, and 
evidence of metals in the photosphere (BBR98; [Bildsten fc Rutledge 2000|) . 



An intensity decrease of 40±8% across 5 yrs in Cen X-4 may well be explained 
as crustal cooling from a long (~100 yr) recurrence time NS ( jKutledge et al. 2001cj ). 
Previous observations have shown shorter term variability (>1 day). In Cen X-4 ( |Campana| 
et al. 199~?D , the intensity decreased by a factor of ~ 3 over 4 daysQin 1997, approximately 



18 years after the previous known outburst. Since, in crustal cooling variability, the 
timescale for variability scales with the time since the most recent outburst, four-day 
variability is much too short a timescale to occur 18 years post-outburst to be explained by 
this mechanism flUshomirsky fc Rutledge 2001 ). 



A second means by which active accretion onto the NS would be indicated during 
quiescence is the presence of metals in the photosphere. For very low accretion rates (an 
order of magnitude estimate is M c <^ClO 13 MQyr" 1 ) a transient NS will have a pure H 
atmosphere because the heavy element settling time is ~ 10 s at the photosphere ( Pildsteri 



|et al. 199~2|) . At accretion rates above this, metals continuously populate the photosphere 
at such densities that their presence may be observable as absorption lines or edges (see 
Fig. @. 

Aql X-l has been detected in X-ray quiescence six times: once with the ROSAT /YUM 
and twice with ROSAT /PSPC QVerbunt et al. 1994| ), once with ASCA ( |Asai et al. lWj ), 



once with BeppoSAX flCampana et al. 1998b|) , and once with Chandra (R01). In all 
but the last (R01), the emission was interpreted as blackbody, although a significant 
power-law component was detected by |Campana et al. ( 1998b )| , which dominated the 



emission above ~2 keV, similar to a component observed from Cen X-4 ( Asai et al. 1996a ; 



Campana et al. 2000|; [Rutledge et al. 200 lc|). R01 examined a number of spectral models, 



discarding them in favor of a H atmosphere spectrum associated with emission from the NS 
photosphere. 

We report here three additional Chandra detections in quiescence. In three of four 



6 Our re-analysis of the data found that the ROSAT/HRI countratc decreased by a factor of 4.1±0.6 over 
4.4 days. 
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Chandra observations of Aql X— 1 in quiescence, we detect neither short term variability, 
nor absorption lines; however, in the fourth observation, we detect both variability and a 
possible absorption line. While other explanations are not excluded, the presence of both 
short-term (<10 4 sec) variability and absorption during the fourth observation, combined 
with the absence of both during the other three observations, is evidence for active accretion 
onto the NS in quiescence. The modest significance (5c) of the absorption line does not 
permit detailed spectroscopic study; however, its presence — should it be confirmed with 
higher S/N observations — opens up the opportunity to measure the redshift at the NS 
photosphere. Simultaneous measurements of the redshift, which is a function of M/R, and 
R would constrain the NS equation of state. 

We begin, in § by reviewing the distance and reddening to Aql X— 1. We then 
describe, in § ^ our observations, timing analysis, and spectral analysis. Section [| discusses 
our interpretation of the spectra and its implications. 



2. Aql X-l: Distance and Reddening 



Callanan, Filippenko & Garcia ( 1999Q showed that the optical counterpart to Aql 
X-l is a faint star 0.46" from the previously mis-identified counterpart. This led to the 
counterpart's identification as a late type star (spectral type K7 to MO) with a quiescent 
magnitude V = 21.6 with reddening of E(B — V) = 0.5±0.1 QChevalier et al. 1999Q . 
For standard conversions^ {A V /E(B - V)=3.1, |Fitzpatrick 19991 ; N n / A v =1.79xl0 21 

Predehl fc Schmitt 1995| ), this corresponds to an X-ray column density of 



cm 



mag 



AT H 22= o.28±0.06 (A H =A H ,22l0 22 cm" 2 ). [I'horstensen et al. (1978)1 note that a nearby (1.4') 
B-type star has an optical reddening E(B — V)=0.73 mag (no uncertainty given), implying 
Ah,22=0.40. These values are comparable to the integrated galactic HI measurements in 
the direction of Aql X-l of A H , 2 2 =0.34±0.01 flDickey fc Lockman 199D| , from W3nHQ; the 
uncertainty was found as cri° / V^A, where cr^ is the dispersion of the N measurements 
within 1 degree, and N = 7). These Ah values are listed in Table |j. These values will be 
used for comparison with the results from spectral fitting in the following section. 

The orbital period has been measured at P or b = 18.95 hr via photometric observations 
both in outburst ( |(Jhevalier fc llovaisky 199"8~5| ; |Garcia et al. 1999] ) and quiescence ( |Welsh| 
et al. 2000]). While Chevalier et al. (|1999|) estimated the distance to the binary as 2.5 kpc, 



7 The formal uncertainty in this conversion is 2%, however there are clear systematic deviations of 25% 
to factors of ~few in some X-ray sources; this may be due in part to the assumed X-ray spectral model. 

8 http:/ /heasarc. gsfc.nasa.gov 
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a re-examination of the arguments producing this distance (which neglected the need for 
Roche-lobe overflow in the primary) finds that it is between 4 and 6.5 kpc (R01); we adopt 
5 kpc as our fiducial distance. This 20% uncertainty is neglected in our quoted uncertainties 
in Roo, which only take into account the statistical uncertainty in our spectral fits. The 
system's orbital inclination is estimated to be > 36° ( [Welsh et al. 20001 ) . 



3. Observations and Analysis 

Aql X— 1 was observed on four occasions, with identical instrumental set-up using 
the Chandra/ ACIS-S detector, backside-illuminated chip. The X-ray source was placed 4' 
off-axis and the read-out limited to 1/8 of the chip area with time resolution of 0.44104 
seconds; these settings mitigate pile-up to < 1% of the detected photons. In all observations, 
source counts were extracted from an 8 pixel radius about the source position; background 
was taken from an annulus of inner and outer radius of 10 and 50 pixels, respectively. 
The spacecraft roll angle differed between each observation. The X-ray source, however, 
fell on the same physical pixels which, due to the spacecraft dither (used to average over 
QE non-uniformity and bad pixels) and larger PSF, spanned ~ 50 detector pixels in both 
directions across the chip. In each of the four observations, the counts were centered near 
(chipx, chipy) = (690,710), with deviations between observations of ±4 pixel; there are no 
known bad pixels within 50 pixels of this location. The expected background countrate in 
the source region was always <1%. Analysis of the first of these observations is described 
elsewhere (R01). An observation list is in Table [|. 

For the analysis we use CIAO 2.2.0.1 with CALDB 2.9; the ACIS-S3 chip calibration 
was improved in CALDB 2.70; analyses of astrophysical sources shows no evidence of 
systematic deviations >15%. However, the effects of a recently reported time dependence 
in the ACIS response is discussed in Sec. 

The first observation triggered as a Target of Opportunity, following the end of a bright 
X-ray outburst (see Fig. P - the brightest observed yet with RXTE/ASM ( |Jain et al. 2000 ; 
Rutledge 2000| ). The other three observations followed after a space of 2, 3 and 4 months 
from the first observation. The countrates in Table |1], in three separate passbands, show 
significant variability between observations 1, 2 and 3; observations 3 and 4 show identical 
countrates in all passbands. In the 0.5-10 keV passband, the countrate decreased by 51±8% 
between observations 1 and 2 (over a period of 81 days), increased by 35±4% between obs. 
2 and 3 (across 32 days), and then remained consistent with constant (±6%) across 28 



http://cxc.harvard.edu/cal/Links/Acis/acis/Cal_prods/matrix/matrix.html 
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days) . 



3.1. Short-Term (<10 4 sec) Variability 

Using the same approach as in R01, we produced power density spectra (PDS) using all 
counts detected (no energy constraints, corresponding roughly to 0.2-10 keV) to search for 
broad-band variability (see [Lewin et al. 1988| for the use of PDS in measuring broad-band 



variability in X-ray sources). We produced Fourier transforms of the time-series data, with 
time resolution of 0.44104 sec ( [Press et al. 1995| ). From these, we produced the PDS (the 



sum of the squares of the Fourier components as a function of frequency), and rebinned 
them logarithmically in frequency. The resulting PDS are shown in Fig. [|. 

Visual examination of the PDS shows no evidence for broad-band power in observations 
1-3 (there is marginally significant variability in observation 3, near 2xl0~ 4 Hz); however, 
there is clear excess power at low frequencies during observation 4. We modeled the data 
as a power-law in frequency, plus a constant (Poisson) level (P{y) = Av~ a + C); the best 
fit model gives a slope a = 1.8±0.4 and an root-mean-square (RMS) variability of 32jJ%, 
integrated between 0.0001-1 Hz (0.2-10 keV). Similar modeling to observations 1-3 (holding 
the power-law slope fixed at the best-fit value of observation 4) provide 3cr upper-limits 
between 18-29% rms during these observations (see Table |2|). For comparison, we show 
the lightcurve of obs. 4 in Fig. |3|. This is the first observation of short timescale (< 10 4 s) 
variability observed from the intrinsic spectrum of a quiescent NS transient (that is, not 
attributable to changes in the intervening X-ray column density). 

We estimate (based on spectral fiting, below) that the power-law component only 
contributes 12% of the counts in the detector in the 0.2-10 keV energy range. This is 
insufficient to produce 32% variability. We conclude that the variability is in the thermal 
component. 

Since the ratio of X-ray to bolometric luminosity log(Lx/£boi)=0.5±0.2 from Aql X— 1 
in quiescence is more than three orders of magnitude greater than observed from coronal 
emission from the analogous RS CVn systems (|Bildsten fc Rutledge 2000]), we exclude the 



possibility that observed variability (at Ri30% of the quiescent luminosity) could be due to 
a coronal flare of the companion, since this would, again, be 3 orders of magnitude greater 
than flares seen from other systems (see also R01). 
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3.2. Spectral Analysis 

Data were extracted, using psextract, into pulse-invariant (PI) spectra. We binned all 
spectral data identically, as shown in their best fit spectral figures (Figs f|a-d); the binning 
was selected so that: (1) low-energy bins are wider than the detector energy resolution; (2) 
the 0.45-0.60 keV bin integrates completely across the O-edge in the detector, over which 
the detection efficiency varies significantly; and (3) the high energy bins have comparable 
signal-to-noise. We analyzed the PI spectra using XSPEC vll.0.1 ( |Arnaud 1996 ). 

Our default spectral model is a H atmosphere spectrum ( |Zavlin et al. 1996| ) plus 
power-law and galactic absorption^]. Other spectral models (blackbody, Raymond-Smith 
plasma, pure power-law, disk blackbody) for this source in quiescence have been examined 
elsewhere (R01), using data from observation 1; note that in the present work, the analyzed 
data extends down to 0.3 keV, while in the previous work, data were not used below 0.5 
keV. 

We first fit the four spectra individually (§ |3.2.1| ), and followed this with combined 
fits of observations 1 and 2 (§ |3.2.2| ), and 3 and 4 (§ |3.2.3| ). This was done to address a 
systematic difference between the two groups of observations. We then attempted (§ |3.2.4| ) 
joint fits of all four observations to determine which spectral parameters changed between 
them. 



3.2.1. Spectral Fits to the Individual Observations 

We began by fitting all four spectra individually with an H atmosphere model 
with galactic absorption and no power-law component. While obs. 1 and 2 are 
statistically acceptably fit (reduced chi-squared statistic xl=2- 05/10 degrees of freedom, 
probability=0.03, xt— 1-10/10 dof, prob=0.36 respectively), obs. 3 and 4 are not 
(^=4.23/10 dof, prob=6xl0~ 6 , ^=6.8/10 dof, prob=10- 10 respectively). Visual 
examination of the spectra from obs. 3 and 4 show systematic excesses in counts above the 
model prediction at energies >3 keV. Such an excess has been observed previously from 
this source (|Campana et al. 1998 b| ) and other qNSs (see § P~2|) , and is typically modeled 



as a power-law component. We therefore include this component in all four spectral fits, 



10 



Neutron star H atmosphere emergent spectral models — and more complicated metallic models — have 
also been calculated by Zavlin et al. (1996), Rajagopal & Romani (1996) and more recently by Gaensicke 
ct al. (2002). The H atmosphere models of these three references agree to within ^few per cent at energies 
relevant to the present work, 0.1-3 keV. 
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to provide upper-limits on its flux in observations 1 and 2. The resulting best-fit spectral 
parameters for the individual observations are given in Table |3|, and figures of the best-fit 
models and data are in Figs. [|a-d. 

There is a systematic difference between the best-fit A% value and Roc values between 
the pairs of observations 1+2 and 3+4. We address this in the following sections, by 
comparing first the spectra of observations 1 and 2, then the spectra of observations 3 and 
4. 



3.2.2. Observations 1 and 2 

As neither of these spectra statistically require a power- law component, we investigate 
if the difference between them is due to A^h, kT e g jQO , or R^. We jointly fit the two 
spectra with an absorbed H atmosphere plus power-law, permitting 1 of the three spectral 
parameters (Ah, /cT e fj j00 , R^) to vary in turn between the two observations. We include 
a power-law component to avoid systematic bias in the derived H atmosphere spectral 
parameters. The resulting best-fit parameters and statistics are given in Table [|. 

The hypothesis of a variable Ah is rejected. A variable R^ is statistically acceptable, 
although not physically motivated in the deep crustal heating model. Finally, the hypothesis 
of a variable fcT e fj j00 is statistically acceptable. We therefore attribute the difference between 
the spectra of these two observations to changes in the kT e g )00 , which decreased from 
125lg 2 eV to IO8I7 1 e V over 81 days. These uncertainties include covariance with other 
spectral parameters (of which Ah and .R^are the strongest). When we hold all other 
parameters at their best-fit values, the uncertainties are 125+1 eV, and 108^2 eV (1 a). 
Note that the best-fit Ah,22=0.43±0.03 is consistent with values found through other 
methods (Table ^). 

Attributing the variability to the magnitude of the power-law component is statistically 
permitted by the data; however, in so doing, the best-fit photon index is required to be 
steep (4.5+0.2) as to necessarily dominate the emission across the full Chandra energy band 
(>95% of the 0.5-10 keV luminosity), the thermal component is not required statistically 
(i?oo<7.2 km at a fiducial value of /cT e g j00 =115 eV), and the value of Ah is a factor of 2 
above the optically implied value. As we know of no physical scenario in which such a steep 
power-law spectrum is produced, and it would require the luminosity to be serendipitously 
consistent with that of a NS with R^^lb km at the observed distance and Ah, we think 
this hypothesis unlikely. We therefore attribute the different intensities of observations 1 
and 2 to a change in fcT eff oo . 
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3.2.3. Observations 3 and 4'- A Possible Counts Deficit in the 0.45-0.6 keV Bin 

We jointly fit observations 3 and 4 with the same absorbed H atmosphere plus 
power-law spectrum. A statistically acceptable fit is found (Table |5|). Thus, Aql X— 1 
observations #3 and #4 are consistent with having the same spectrum and intensity. The 
best fit Nft t 22 (0.60 ±0.05) is greater than that found by other methods (Table f|), including 
the value found from joint spectral fitting from observations 1 and 2. The value of Roo is 
marginally larger than in observations 1 and 2 (46t^ km vs. 17.2l|; 6 km). In addition, 
the value of Roo is also larger than that found in a previous analysis of observation 1 alone 
(R01), in which we had used counts with energies >0.5 keV, and that found in analyses 
of ROSAT and ASCA data in which N-r was held fixed at the optically implied value 
( IRutledge et al. 1999|) . 

When we hold N^22 and R^ fixed at the best-fit value from observations 1 and 2, the 
best-fit spectrum is marginally acceptable (x^=l. 71/25 dof, prob=0.015). Examination of 
the spectral model and data (Fig. |5|) indicates a deficit in the 0.45-0.6 keV energy bin: 
the number of detected counts appears significantly below that predicted by the best fit 
model. When the 0.45-0.6 keV energy bin is removed from the fit, the new best-fit model 
is significantly improved (x^=0. 76/23 dof; prob=0.78). This suggests a deficit of counts in 
the 0.45-0.6 keV range in observations 3 and 4 (largely, the latter). 

The best-fit model (ignoring this energy bin, with A^h,22 and R^ held fixed at the 
best-fit observation 1 and 2 values) predicts 20 and 25 counts in this bin for obs. 3 and 4 
respectively (the difference in counts is due to different integration times). Only 14 and 10 
counts are observed in observations 3 and 4. For an average countrate of 25 counts/bin, the 
Poisson probability of detecting <10 counts is 0.06%. Combining the two observations, the 
probability of detecting <24 counts (with an average of 45 counts/bin) is 0.04%. No such 
deficit is apparent in the spectra of observations 1 or 2. 

3.2.4. Joint Fitting 

To examine which spectral parameters are different between the four observations, we 
jointly fit all four spectra simultaneously. Because of the similarity between observations 
3 and 4, we treat their spectra as being identical (i.e., all parameters for observations 3 
and 4 are treated as the same), so we are examining the difference among three spectra, 
observations 1, 2 and 3+4. Table ^| contains a summary of our findings. 

We fit the data with the absorbed H atmosphere + power-law spectrum and limited 
our fits to the 0.6-9.0 keV energy range, to avoid affects associated with the deficit below 
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0.6 keV in observation 4. As may be surmised from the different countrates, using a single 
spectrum with all five parameters (A Hi 22) Roc, a , Fx, P u an d &^eff,oo) f ree but the same for 
all 3 spectra, the best-fit is unacceptable (x^=5. 52/43 dof; prob=10~ 28 ). We then allowed 
each of the spectral parameters to vary between observations. 

1. A fit with kT e g }OQ allowed to vary is statistically acceptable, for values of 
fcT effj00 =121^ 3 eV, 105±|? eV, and HOif eV for observations 1, 2, and 3+4 
respectively. When we hold all other parameters fixed at their best-fit values (to 
exclude covariance in the uncertainties) the values are 121±1, 105±1, and 110±1 eV. 

2. When we permit only Ah to vary between the three spectra, the best-fit is 
unacceptable (x^=3. 05/41 dof; prob=2x 10~ 10 ). Thus the changes between the 
spectra cannot be due only to a changing A H . 

3. Although not physically motivated, when we permit only to vary between the 
three spectra, we obtain an acceptable fit (x£=l. 13/41 dof; prob=0.26). The range of 
acceptable Roo values is small: 12.5-17.6 km. 

4. We then fit the data with only Fx, P i allowed to vary. The best fit was not statistically 
acceptable (Xy=l-60/41 dof; prob=0.0085). When we permit both a and F x>p i to 
vary, the fit becomes acceptable; the changes between the power-law slopes are large, 
however, and the best-fit R^ is larger than expected from theory, and Ah is larger 
than that found through radio and optical observations. It is unlikely, therefore, that 
a variable power-law component is responsible for the different spectra. 

The spectra of observations 1 and 2 are acceptably fit without a power-law component, 
and the combination of the three spectra are acceptably fit with only a changing kT e g jOQ . 
This motivated a fit with the following parameters: no power-law component in observation 
1 and 2; A H held fixed at the radio measured value (A Hi 22 =0.34); R^ held fixed at 13.0 km 
(a value which nearly all theoretical EOSs can produce; [Lattimer fc Prakash 200l| ); and 
kT e s t00 varying. The best-fit is acceptable (Xy=l-30/43 dof; prob=0.11). When we now 
expand the energy range of this fit to include 0.3-0.6 keV, however, the best fit becomes 
unacceptable (xl=2. 63/51 dof; prob=2x 10 -9 ). Examination of the spectrum shows that 
the observed countrates are systematically lower than the predicted countrates in the 
0.3-0.6 keV range. When we let both Ah and R^ float (with the same value in all three 
spectra), the best-fit spectral model is still not statistically acceptable (xl=l. 62/49 dof; 
prob=0.0038). This unacceptable spectral model is shown in Fig. |6|. The most discrepant 
point is the 0.45-0.6 keV bin from observation 4. When this point is removed, a new best-fit 
model is acceptable (x^=l-26/48 dof; prob=0.105). 
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Since the poor fit to the full data set can be attributed to the single low-energy bin in 
observation 4, one possibility is that N-r changed between observations 1+2 and observation 
3+4. If we permit the values of Nn and kT e g ]OC to change between observation 1, 2 and 
3+4, however, the spectrum is still not acceptably fit (x^=l-66/47 dof; prob=3x 1CT 3 ); 
therefore, the deficit is not explained by a changing JVh. 

In summary, the differences in the three spectra are acceptably explained as entirely 
due to a changing fcT e fj j00 of the surface thermal emission. The differences cannot be 
explained as entirely due to either a changing power-law flux or to a variable Nr. While a 
combination of variability in a and Fx, P i provides a statistically acceptable fit, the required 
variation in the power-law slope a is large and the resultant best-fit values of Roo and 
]Vh are larger than we would a priori expect. We therefore conclude that the simplest 
explanation for the variation in the spectrum of Aql X— 1 in quiescence is due to a change 
in /cT e ff j00 of a pure H atmosphere. 



3.3. Re-examination of the Counts Deficit in Observation 4 

We now turn our attention to the apparent deficit in counts in the 0.45-0.6 keV bin 
of observation 4. In particular, we examine whether it is possible to explain the deficit by 
some instrumental or analysis effect that is unrelated to the NS emission. 

When we exchange the response matrix from obs. 4 to that from obs. 1, there is no 
effect on the spectral fit (the focal plane temperatures were the same in all observations, 
and the X-ray source was focussed at the same location on the S3 chip). The deficit is 
therefore not caused by a change in the calibrated spectral response particular to obs. 4. 

What if the effective area of the 0.45-0.6 keV bin is simply different from that 
represented in CALDB v2.9? Could we simply change the expected number of counts 
(on the basis of the best-fit spectrum and the calibrated response) by the same constant 
factor, which would make all observed 0.45-0.6 keV countrate consistent with the expected 
number? The calibrated response is based on detailed knowledge of the telescope+detector 
system. We took the response matrix to be completely diagonal, such that a photon of 
energy E n would be detected as a count with energy E n . While the system response is 
considerably more complicated, with significant non-diagonal elements, this zeroth order 
approximation can reveal if change in the 0.45-0.6 keV effective area can account for the 
observed discrepancy, or if higher order changes would be necessary. 

In Table we list the number of counts detected in the 0.45-0.6 keV energy bin 
during each observation i, and the number of counts predicted to be in the bin (/ij) by our 
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best-fit spectral model (in which the 0.45-0.6 keV bin was ignored in observation 4). We 
then calculated a probability for the factor (/) by which we change the diagonal response, 

Prob(f) = nf =1 P Poisson (/ x Hi,Xi) 

where -Pp isson(/W, x) is the Poisson probability of observing X{ counts when fa is the 
expected average realization. We find, for our observed and predicted values of Xi 
and fa, that Prob(f) is maximal when / = 0.75 (that is, the observed counts would 
be most consistent with an effective area in the 0.45-0.6 keV bin which is 75% of the 
telescope+detector effective area represented in CALDBv2.9). We performed a Monte 
Carlo simulation, with 4 Poisson deviated counts, with mean counts of f x fa, calculating 
ProfeMonte-Cario(/ = 0.75) for 10 6 such realizations. We found that we would produce a 
value of -Pro6MonteCario(/ = 0.75) < Prob(f = 0.75) 2.0% of the time (that is, even the most 
likely value of / would produce the distribution of observed only 2% of the time). We 
thus marginally exclude the counts deficit as due to a different effective area in the 0.45-0.6 
keV bin from that in CALDB 2.9. To attribute the observed discrepancy to a response 
inaccuracy, there must either be a significant difference between a non-diagonal element as 
modelled and its true value, or the diagonal element response in the 0.45-0.6 keV bin (and 
only that bin) may be time dependent. 

We also examined the spectrum in PHA space as well as PI space, and find the same 
deficit of counts in observation 4. A statistically anomalous background subtraction (i.e, 
more than an average number of background counts in this spectrum) is not responsible for 
the deficit either; when we neglect the background entirely (which we estimate to be ~9 
counts of the 1137 counts observed in this spectrum), the deficit remains. 

We therefore marginally exclude the deficit of counts to response matrix uncertainty 
under a diagonal assumption, and exclude either anomalous background, or binning in 
PI space vs. PHA space. The deficit could be due to non-diagonal elements in the 
redistribution matrix being different than those in CALDB v2.9, or to an unknown 
time-dependency in the diagonal elements of the redistribution matrix. 

3.4. Time Dependent Response Below 0.7 keV in ACIS-S3 

When this work was largely complete, it was reported that there is a time-dependency 
in the response of the ACIS chips below 0.7 keV 0, which appears as a decrease in effective 



http:/ /asc. harvard.edu/cal/Links/Acis/acis/Cal_projects/index. html 
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area, the magnitude of which has increased with time since launch. The magnitude, exact 
time dependency, and energy dependency are not fully known at present, and it is beyond 
the scope of the present work to provide a full time- and energy- dependent re-calibration 
effort. However, the discussed possible causes (deposition of contaminants on the CCDs or 
the optical blocking filter) suggest the effect would be monotonicaly increasing with time, 
perhaps at a constant rate. 

We examined briefly the magnitude of this effect, to determine its possible effects on 
our conclusions. We used two observations of the spectrally soft super-nova remnant SNR 
E0102-72.3, taken 2000 May 28 (ObsID 141, ~6 months prior to Obs. 1 of Aql X-l) and 
2001 Dec 6 (ObsID 2844; ~8 months after Obs. 4 of Aql X-l). These were taken with the 
same focal plane temperature (— 120F) but at different locations on the ACIS-S3 chip (141 
was 2.2' off-axis, while 2844 was only 1' off axis). We extracted data within 39" of the SNR 
center, and produced pulse-invarint (PI) spectra with resolution (10-15 eV) greater than 
the intrinsic energy resolution of ACIS-S3 between 0.3-1.1 keV (120-140 eV). We modelled 
the first observation with a model meant to be parametric (wabs * (c6vmekl + gauss)), 
which with a best fit of x^=3.30 for 32 dof was not statistically acceptable; however the 
most deviant PI bins were not systematically offset from the model (that is, the scatter 
appeared random as a function of energy), and were at most 3a from the model. This 
is adequate for our purpose of investigating systematic response differences at the energy 
resolution of the detector. 

We compared the second spectrum obtained 557 days later with the best-fit model 
of the first spectrum (see Fig. ||). As can be seen from the figure, the response appears 
to be diminished by ~50% between 0.3-0.5 keV, with the magnitude of the discrepancy 
decreasing between 0.5-0.7 keV. Thus, from our comparison, we see an apparent change 
in the detector response, which may be modelled as an energy-dependent decrease in the 
effective area below 0.7 keV. 

We investigate the effect of this on our conclusions regarding the broad-band spectrum, 
by using data only in the energy range 1-10 keV. First, as can be seen from Table [I], 
the 1.0-2.5 and 2.5-8 keV countrates vary between observations, and so cannot be 
explained as due to the <1 keV response. Joint spectral fits using all four observations 
still find that the variation cannot be explained as due to a change in iVn or power-law 
component normalization alone (both with values of xl with corresponding probabilities 
of <10~ 3 ). Since we are simply removing data from consideration, all those combinations 
of parameters which produced acceptable spectral fits using 0.3-10 keV still produce 
acceptable fits using 1-10 keV, but with larger error bars. Thus, the conclusion that the 
observation-to-observation variability is explained by a change in kT e s tOQ remains. We 



-14- 



expect the time-dependent calibration of the sub-1 keV range, when it is fully known, will 
change the resulting spectral parameters we have measured here, but not our qualitative 
conclusions. 

Finally, as regards the apparent counts deficit, the change in effective area does affect 
the energy range of interest (0.45-0.6 keV). We performed the following check. First, using 
the E0102-72.3 observations, we included a multiplicative spline model in XSPEC to 
account for the change in calibration; this was done by first finding the best-fit parametric 
model to the May '00 observation only, and using the same parameters held frozen for 
the Dec '01 data, but including a spline model, with Estart=0.3 keV, Eend=0.7 keV, 
Yend=1.0, YPend=0.0, and the values of Ystart and YPstart permitted to vary. The 
best-fit spectrum was unacceptable (x^=15.3, 52 dof), although the gross change in response 
was well accounted for, with best-fit values of Ystart=1.5 and YPstart=— 20. We then 
refit Obs. 1-4 with the spectral model which had been previously unacceptable due to the 
counts deficit (kT c s tOQ varies between observations, Ah and Roo permitted to float, and a 
power-law spectral component for Obs. 3 and 4 only), but including the spline component 
for Obs 3 and 4. This models the change in response as essentially instantaneous between 
observations 2 and 3. The best-fit joint spectrum is statistically acceptable, and the deficit 
is no longer significant (x^=1.23 / 49 dof; prob=0.13). The accuracy of this correction, 
however, depends in detail on the time- dependence of the change in response, and whether 
or not it can be modelled as a change in effective area as a function of energy. This requires 
much more detailed modelling of the response than we can do here. We conclude, however, 
that the deficit of counts can be due entirely to the change in detector response; this requires 
confirmation through detailed study and modelling on the time- and energy- dependence of 
the sub-0.7 keV detector response. 

4. Discussion and Conclusions 

Following an outburst of Aql X— 1 during November 2000, we took four Chandra/ ACIS-S 
snapshots over a span of 5 months. With this series of spectra, we find the following. 

1. All four spectra are acceptably fit with H photosphere models, with radii consistent 
with that of a NS, at the distance and reddening of Aql X-l. In two of the four 
spectra, an additional power-law component is required to acceptably fit the spectrum 
at high energies (>3 keV), as has been found previously for this and other qNSs. 

2. The intensity decreased by 50% over the first three months following the outburst. 
It then increased by 35% over the next month and then remained constant for the 
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last month spanned by our observations. The change in the spectrum between the 
observations cannot be explained exclusively as a change in the power-law component, 
nor in the value of Nr. A change in the photospheric effective temperature is required 
to account for the change. 

3. No short timescale (< 10 4 s) variability was found in the first three observations; 
the fourth had significant variability (32tg% rms, 0.2-10 keV) on timescales 10 3 -10 4 
sec. As the power-law component contributes a smaller fraction (~12%) to the 
total counts, this magnitude in variability cannot be attributed to the power-law 
component. 

4. A deficit of counts in the 0.45-0.6 keV band was noted during observation 4. This 
deficit was not present in the other three observations. Follow-up spectroscopy and 
detailed study of the change in the detector response is required to determine if the 
deficit is intrinsic to the source. 

5. A best-fit spectral model, excluding the 0.45-0.6 keV bin in observation 4, gives a 
value of AH,22=0.42tQ Q3 - consistent with the optical and radio implied value; the best 
fit i?oo=15.9t2]9 km (neglecting the 20% uncertainty in the 5 kpc distance), consistent 
with the radius of a NS for most proposed NS EOS ( [Lattimer fe Prakash 200l"| ). 



The increase in temperature between observations 2 to 3, if significant, is interesting. 
While accretion can clearly affect the thermal spectrum (e.g., Zampieri et al. 1995), 
it is also possible that differential sedimentation of ions in the atmosphere can change 
the surface effective temperature by modifying the opacity at intermediate densities, 



10 g cm <z p <^ 10 g cm , ( [Brown et al. 2002|) . Further detailed calculations are 



required to determine the amplitude and timescale of this effect. 

At present, there is no known mechanism associated with crustal heating that can 
produce the short-timescale (<10 4 sec) variability we observe here. The minimum timescale 
for post-outburst intensity variability in quiescence from deep crustal heating is days 
to months (the thermal diffusion timescale from the depth of heat deposit), and should 
scale roughly as 5t, where St is the time since a (short) outburst ends ( |Ushomirsky fc 



Rutledge 200 1|) . For NSs with "normal" cores, the magnitude of variability is a few per 



cent; if more exotic neutrino cooling takes place in the core, however, the magnitude of 
variability in quiescence can be a factor of ten or more. In the absence of variability due 
to post-outburst crustal cooling, the only other timescale for variability is the timescale of 
core cooling, which is typically 10 5 -10 6 yr. In contrast, intensity variability due to accretion 
onto the compact object can in principle take place on the dynamical timescale, seconds 
or shorter. However, while the observation is in qualitative agreement with accretion onto 
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the compact object, it remains a major puzzle why the magnitude of luminosity should 
be quantitatively comparable to the luminosity predicted previously from deep crustal 
heating R01. The absence of a quantitative prediction for the magnitude of accretion onto 
the compact object in quiescence further complicates comparison between this observation 
and an accretion theory. For the present, we conclude that the intensity variability on 
<10 4 sec timescale can be explained by accretion onto the compact object, but the implied 
coincidence of the accretion luminosity being comparable to the predicted deep crustal 
heating luminosity, as found by R01, remains a challenge to theory. 



4.1. Uncertainties in the Measurement of R, and Prospects for a Useful 

Constraint on the EOS 

Throughout this paper, we use the word "radius" to mean the measured value i?^, the 
effective emission area size observed from an infinite distance. In practice, what we measure 
using spectral analysis is the angular size of the source R^o/D, where D is the distance to 
the NS. We assume D =5kpc, a value which is at present uncertain by ±20% (R01); we 
neglect this uncertainty throughout the present work, addressing only the statistical and 
spectral uncertainty, which can be greatly improved upon with higher signal-to-noise X-ray 
spectra. The distance to Aql X— 1 can be accurately measured in the future using the Space 
Interferometry Mission (SIM; |Unwin 

199D, which can achieve ±2% accuracy at 5 kpc for an 
object with V < 17m (Aql X-l obtains V = 14.8 in outburst ; [Liu et al. 2001|) . Such precise 



measurement will effectively eliminate distance as an uncertainty in R^. Indeed, using 
spectral analyses of qNSs in globular clusters ( |Grindlay et al. 2001| ; [Rutledge et al. 2001b 



jin't Zand et al. 2001|) , the distances to which can be measured with ±2% precision using 
Hipparcos dwarfs ( |Carretta et al. 2000 ), is one approach being pursused to obtain accurate 



Roo measurements prior to the launch of SIM. 

We have assumed throughout this work that i?oo does not change. In the deep crustal 
heating scenario (BBR98), luminosity originates from the NS core, and the NS star has 
a kT e g yOQ which is isotropic at the photosphere, which we assume throughout this paper. 
This scenario, however, is not entirely consistent with the short-term (seconds-months) 
variability which is observed. If the photosphere were to have temperature anisotropies - for 
example, which may occur if the NS were powered exclusively by non-spherically symmetric 
accretion - this would change the physical meaning of Roo, which would then dependent on 
the details of the photospheric temperature distribution, as described for accreting NSs at 
much higher luminosities than are relevant here QPsaltis et al. 20001) . 



The value of R^ is related to the physical radius R by R OQ =R/ Jl — 2GM/(Rc 2 ) 
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where G is the graviational constant, M is the NS mass, and c is the speed of light. It is 
the value R which has recently been shown to place a useful constraint on the NS EOS 
( [Lattimer fc Prakash 20(JT| ), if R is measured to <10% accuracy. The value R cannot be 
derived from i?^ in the absence of a known NS mass M. If one assumes a range of values 
between O.8-2.8M , this introduces an additional systematic uncertainty of ~50% in R, 
which is greater than that needed to usefully constrain the EOS in the prescription of 



Lattimer fc Prakash (200"T]| . Thus, precise measurements of R^ will not usefully constrain 
the NS EOS in the absence of a measured NS mass or photospheric redshift; however, if a 
phenomenological relationship between R^ and the pressure at nuclear density for a range 
of EOSs can be found, as found for R by [Lattimer fc Prakash (200T]| , it may be that Roo 
will usefully constrain the NS EOS. This possibility is presently under examination. 



4.2. The Power Law Component 

A power-law spectral component is required to fit the spectrum above 3 keV in 2 of the 
four observations, and it is consistent with being the same magnitude and slope in all four 
observations. 

A high energy power-law tail has now been observed to dominate the X-ray spectrum 
above 2-3 keV in five qNSs: Aql X-l, Cen X-4 QAsai et al. 1996U |Campana et al. 200q ; 



Rutlcdge et al. 200q) , KS 1731-260 flWijnands et al. 2001| ; [Rutledge et al. 2002|) , NGC 



6440 flin't Zand et al. 200 1|) and 4U 2129+47 ( Nowak 20021) . In all cases, this component 



comprises ~10-40% of the 0.5-10 keV luminosity. Importantly, in all those power-law 
components which have been detected, the slope appears to be either inverted or flat. A 
cutoff at higher photon energies (above 3 keV) is required for the luminosity not to diverge. 

It has previously been suggested that the component may be due to accretion onto 
the NS magnetosphere (|Campana et al. 1998a| ); however, no detailed a priori model 



(including an energy budget, estimation of the accretion rate onto the compact object 
or magnetosphere) exists. If the power-law emission were powered by accretion onto the 
magnetosphere, however, it would be coincidental that five different sources (each having 
different B fields, spin rates, and perhaps accretion rates onto the magnetosphere) would 
have comparable luminosity ratios of the power-law component to the thermal component, 
which originates from the NS surface. The same coincidence would apply to non-thermal 
emission from pulsar shocks as well QTavani 199 1| ; (3tella et al. 1994] ; [Menou et al. 199"9| ). 



One possible suggestion is that both the thermal and power-law spectral components 
are powered by low-level accretion. Using the assumption of spherical flow onto an 
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unmagnetized, nonrotating neutron star, and assuming that only Coulomb forces acted 
to decelerate the infalling protons, ^ampieri et al. (1995)| showed that no similarly hard 
spectral component was expected at accretion rates comparable to that required to power 
the quiescent thermal emission. Ueulel et al. (2U0T)| similarly found that, for accretion via 
an advective flow with virialized protons, the emergent spectrum is slightly harder than 
thermal, though the hardened component still decreases in vF v . It is therefore unclear that 
the emergent spectrum of a low M NS would produce a hard tail such as observed here. 



4.3. A Deficit of Counts in the 0.45-0.6 keV Bin in Observations 3 and 4 

We found that, compared with a spectral model interpolation in the 0.45-0.6 keV 
energy range, there is a deficit of counts significant at the 99.96% level, during observations 
3 and 4. Such a deficit is not present in observations 1 or 2. There are at least a few 
possible explanations for this. 

It may simply be that the deficit is a statistical anomaly, which will not be confirmed 
with deeper spectroscopy. We may be able to address this possibility with our approved 
30ksec AO-1 GO observation of Aql X-l with XMM-Newton. 

A simple decrease in the effective area in the 0.45-0.6 keV energy bin, by a factor 
of 0.75, is marginally excluded (2% probability) as the source of the counts deficit. The 
deficit may be due, however, to non-diagonal elements of the photon energy redistribution 
matrix which are different from those in the present calibration model. Moreover, a recently 
reported time-dependence in the ACIS response below 0.7 keV can expalain the entire 
deficit, when the time dependence of the change in response is modelled as having taken 
place between Feb and April 2001 (observations 2 and 4). Until this response is properly 
modelled, it cannot be neither fully implicated nor exonerated as being responsible for the 
counts deficit. Thus, an unaccounted-for discrepancy in the detector response remains a 
possible explanation. 

Finally, it may be that the deficit of counts is due to the presence of metals. This 
would be consistent with the observational scenario expected to produce such features 
( [Brown et al. 1998 ), in which active accretion onto the NS surface replenishes the surface 



metallicity. The detected intensity variability (32tg% rms) argues for active accretion 
during this observation as well (such variability may also be present during observation 3); 
and the luminosity of Aql X— 1 during these observations is consistent with this scenario as 
well. 



We show in Fig. [7] the emergent spectra of a solar-metallicity atmosphere from a NS 
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with the effective temperature of Aql X— 1, rotating at 550 Hz ( Zhang et al. 1998| ), at three 



different viewing angles with respect to the rotation axis (the orbital plane of Aql X— 1, 
is a > 36°). The absorption features due to metallic absorption lines are apparent in the 
pole-on (a = 0°) spectral model. The most prominent feature, smeared by NS rotation, 
is near 0.9 keV, which is a mixture of OVII absorption, plus other metal absorption lines. 
In a rough approximation: to redshift this feature to 0.45-0.6 keV requires a redshift 



l + z = g -i = i/^i _ 2GM/(Rc 2 )=l. 5-2.0, which, for a NS of R O0 =Rg~ 1 =15.9 km, implies 
a coordinate radius R = 8 — 10.5 km, and a mass M ^2.0M Q . This is a mass/radius 
estimate of scale, which shows consistency with theoretical expectations, but is not a precise 
estimate due to the large uncertainty in the redshift factor, in addition to the systematic 
uncertainty in the correct identification and reality of the feature in question. 

If future observations confirm the presence of this absorption line in Aql X— 1 in 
quiescence, it can be exploited. By measuring the redshift of such lines (g" 1 ) simultaneously 
with the emission area from the spectrum (oc R 2 g~ 2 ) the NS mass and radius can be 
determined, providing a strong constraint on the equation of state of dense matter ([Lattimer 
fe Prakash 2001) . 
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Fig. 1.— Lightcurve of Aql X-l, outburst through quiescence. RXTE / ASM points 
(squares) used the ASM countrates, assuming 1 ASM c/s=4.4x 1CT 10 erg cm -2 s" 1 (0.5- 
10 keV; appropriate for a kT = 5 keV thermal bremsstrahlung spectrum corrected for 
absorption), and the ACIS-S points (triangles) assume the spectrum for observations 3 and 4 
(1 ACIS-S c/s=1.0xl0- n erg cm" 2 s- 1 , 0.5-10 keV). Only ASM detections of > Aa are used. 
The uncertainty in luminosity due to the countrate is shown; the uncertainty in luminosity 
does not include the systematic uncertainty in the bolometric correction. 

Fig. 2.— PDS of observations 1-4. Points are data, normalized according to Leahy 
et al.(1983). Error bars and upper-limits are la. Solid lines are the best-fit broad-band 
models (see text). RMS variability model values are 0.0001-1 Hz, for all counts data 
detected (roughly, 0.2-10 keV). Observations 1-3 give only upper-limits for variability, while 
in observation there is a clear detection of variability, with 32^6% rms. 



Fig. 3. — Lightcurve of observation 4. 



Fig. 4. — Best fit spectra for observations 1-4 (parameters are in Table ^]). Crosses are 
the observed data. The dashed line is the thermal component, the dashed-dotted line is the 
power-law component, and the solid line is the summed spectral model; the spectral models 
have been corrected for galactic absorption, (a) Obs. 1. (b) Obs. 2. (c) Obs. 3. (d) Obs. 
4. 

Fig. 5. — (Top Panel) The observed spectra for observations 3 (open squares) and 4 (open 
circles), and the best-fit model in which the column density is held fixed at A^h,22=0.38 - 
the best-fit value from observations 1 and 2, as well as comparable to the optically implied 
value (iVH,22=0.35). The spectrum is statistically unacceptable, due to the deficit of counts 
in the 0.45-0.6 keV bin. (Bottom Panel) The residuals \= (data-model) /uncertainty. 

Fig. 6.— (Top Panel) Observations 1-4, with an (unacceptable) best spectral fit (see 
text), due to the low countrate in the 0.45-0.6 keV bin of observation 4. Observation 
1 (stars), 2 (triangles), 3 (squares), and 4 (circles). (Bottom Panel) Residuals x= (data- 
model) / uncertainty. 

Fig. 7. — The emergent (top of the atmosphere - that is, un-redshifted) X-ray spectrum for 
a NS atmosphere (1.4M , 10km) of solar metallicity including the effects of a NS rotating 
at 550 Hz, as observed from Aql X— 1, at three viewing angles (a) relative to the direction 
of NS rotation. The solid line, in which absorption features are apparent, is for an angle 
a = of the NS rotation axis relative to the line of sight. The broken line is for a = 30°, 
and the smooth solid curve us for a = 90°. Z indicates the percentage of metals, assuming 
the mixture of Grevesse and Noels (1993). For Aql X— 1 the best estimate of the rotational 
inclination is a > 36°. 

Fig. 8. — Top Panel: Comparison of two X-ray spectra for the soft SNR E0102-72.3 between 
0.3 and 1.1 keV, taken 5 May 2000 (crosses) and 6 Dec 2001 (boxes). A parametric model 
(solid line) is fit to the the May '00 observation. The Dec '01 spectrum is systematically 
below the May '00 at energies 0.3-0.7 keV. (Bottom panel) Ratio of the best-fit May '00 
parametric spectral model to the May '00 (crosses) and Dec '01 (boxes) data. The systematic 
difference in spectra shows a strong energy dependence: approximately 50% near 0.3-0.5 keV, 
increasing to consistency (100%) near 0.7 kev. Comparison of the increases and decreases in 
the PI spectra also indicate a possible gain shift between the two observations, by ~10 eV. 



"i — i — i — i 1 — i — i — i — i 1 — i — i — i — i 1 — i — i — i 1 — i — i — r 



37 



o 

CD 

w 
u 

CD 



§ 35 



3 

g> 34 



33 



J I I I I I I I I I I L 



J I I I I L 



50 100 150 

Time (Days since 2000 Nov 1.0) 



Figure [l| 



TTTTTTl 1 I I llllll 1 I I I Hill 1 I I MINI TTT 



TTTTTT1 1 I I llllll TTT 



| 1 I I lllllj TTTTI 

Obs. 2 
rms<26% 



10 



I 



Obs. 1 
rms<18% 



* * ■ 



10 - 




Obs .3 
rms<29% 




Obs. 4 
rms = 32!|% 




1 + t * ' 



I I I I Ullllll I I I I I UN* I I I I I 

10" 4 10" 3 0.01 0.1 1 10" 4 10" 3 0.01 0.1 1 



Freq. (Hz) 



Figure |I| 



0.16 



Obs. 4 



0.14 - 



\0.12 
o 



0.1 - 



0.08 



1 



1 



1 



1 







2000 4000 6000 8000 

Time (sec since start of Obs. 4) 



10 4 



Figure || 



•d 
a 

CO 

o 



> 

o 

05 
I 

CO 
I— I 

CO 



CQ 
I 

CQ 
i— i 

O 

< 

as 
U 

-d 

a 



I 

X 




ro 



TOO 



8 



_0T t _0T Z 



2-t-9- 
X 

Figure || 




A931/09S /S^UTIOO 



Figure 




1 

Energy (keV) 



Figure [j] 



i i i — r 




+ 



■ ■ ■ L 



1 1 I 1 1 1 1 I 






? 



+ UL 



s 



+ t 



J L 



I .... I 



08 OT 9 2 T 90 

A93f/0aS /S^UTIOO 



T 90 

Figure 



Table 1. Log of Chandra / ACIS-S Observations 



Obs # 


Start Time 
(TT) 


Exposure 
(sec) 


c/ksec 
(0.5-10 keV) 


c/ksec 
(0.3-1.0 keV) 


c /ksec 
(1.0-2.5 keV) 


c/ksec 
(2.5-8 keV) 


Orbital Phase" 

0orb 


1 


2000-11-28 10:52:39 


6628 


183±5 


56±3 


120±4 


8.3±1.1 


0.02-0.15 (±0.01) 


2 


2001-02-19 11:25:52 


7787 


94±3 


32±2 


60±3 


3.9±0.7 


0.20-0.36 (±0.02) 


3 


2001-03-23 20:09:11 


7390 


127±4 


39±3 


81±3 


8.2±1.1 


0.19-0.34 (±0.02) 


4 


2001-04-20 11:53:32 


9245 


123±4 


38±2 


76±3 


9.9±1.0 


0.22-0.39 (±0.02) 



Note. — a Orbit al phase relative to minimum light (inferior conjunction of the secondary), ephemeris from 
orbital period from Chevalier & Ilovaisky (1998a) 



Garcia & Callanan (1999) 



Table 2. Intensity Variability 



Obs. # 


a: 


Total Counts 


RMS (%) 


1 


(1.8) 


1237 


<18 


2 


(1.8) 


755 


<26 


3 


(1.8) 


956 


<29 


4 


1.8±0.4 


1161 


32±1 



Note. Values in parenthesis are 

fixed. RMS is percentage root mean square 
variability, integrated between 0.0001-1 Hz, 
in excess of the Poisson level; error bars are 
la, upper-limits are 3cr. 



Table 3. Chandra Best-Fit Spectral Parameters for Individual Observations 



Obs. # 


-?Vh,22 


(eV) 


■Roo 

(km/(D/ 5 kpc)) 


q: 





(%) 


Flux" 
(0.5-10 keV) 


xl (prob) 


1 
2 
3 
4 


0.46±0.04 
0.37±0.06 

u - 08 -o.io 
66 +0 - 04 


118± 8 6 
120±20 
S3±\l 
6T±S 


20+ 3 
13±f 
44±17 


(1.0) 
(1.0) 
1.2±0.6 
1.9±0.5 


<2xl0- 5 
<1.2xl0~ 5 
3^2 xl0~ 5 
lto.5 xlo ~ 5 


<17% 
<22% 
18% 
24% 


18 

8 

19 
23 


0.83/10 (0.07) 
0.24/9 (0.99) 
0.98/9 (0.45) 
1.3/9 (0.25) 



Note. — a X-ray fluxes are corrected for galactic absorption, in units of 10 — 13 ergcm~ 2 s -1 (0.5-10 keV). N„i is the 
power-law component normalization in phot keV - 1 cm" 2 s" 1 at 1 keV. f p i is the fraction of the the Total Model Flux which 
is accounted for by the power-law component (these are typically uncertain by 50-100%). Upper-limits are 90% confidence, 
uncertainties are la. Values in parenthesis arc held fixed. Assumed source distance D=5 kpc. 



Table 4. iYn in the direction of Aql X— 1 



Ref. 



Method 



Value 



Callanan et al. (1999)| 



Dickey fc Lockman (19901 (W3NH) 
TTiorstensen et al. (1978)| 



Obs. l+2 a (0.3-9 keV) 
Obs. 3+4 a (0.3-9 keV) 



Optical Photometry 0.29±0.06 

Weighted average integrated 21 cm emission 0.34±0.01 

Optical Photometry of nearby B star 0.40 

X-ray spectral modelling 0.43±0.03 

X-ray spectral modelling 0.61±0.06 



Note. — The value for observations 3+4, without accounting for the apparent 0.45-0.6 keV deficit, 
appears discrepant 



Table 5. Spectral Parameters for Observations 1+2 and 3+4 









Roo 




N p i f p i 




Obs. # 


^H,22 


(eV) 


(km/(D/ 5 kpc)) 


Q 


(phot/cm-^- 1 1 keV) (%) 


X? (prob) 



Observations 1+2; TVjj varies 

1 n/a n/a n/a n/a n/a n/a 5.0/22 (lxlO -13 ) 

2 n/a " " n/a n/a n/a " 
Observations 1+2; R 

oo varies 

1 0.41+0.03 122tg 1 18.0+3.5 0.4+0.9 <4.5xl0" 5 <11 0.83/22 (0.69) 

2 " " 12.9+2.6 " " <21 
Observations 1+2; kT c g ^ varies 

1 0.43+0.03 125lg 2 17.2l|;g 0.9+1.0 <9xl0" 5 <21 0.67/22 (0.87) 

2 " 108+" " " " <35 
Observations 1+2; Fx pi varies 

1 0.83+0.04 (115) i?oo<7.1 4.5+0.2 (2.1±0.3)xl0" 3 >97 1.2/22 (0.36) 

2 " " " " (1.1+0.2) XlO" 3 >95 

Observations 3+4 

3+4 0.60+0.05 81+12 46±^ 1.5+0.5 (5+2) xlO" 5 21 1.2/23 (0.25) 

3+4 (0.43) 113+1 (17.2) 0.8+0.3 (1.7tJ;°)xl0- 5 26 1.71/25 (0.015) 

3+4 6 (0.43) 114+1 (17.2) 0.7+0.4 (1.5+H)xlO _5 26 0.76/23 (0.78) 



Note. — Uncertainties are lcr, upper-limits are 90% confidence. Assumed distance is 5 kpc( Rutledge et al. 2001a). Notation 
"n/a" (not applicable) is used where spectral parameters were used in a fit which is not statistically acceptable. Upper-limits 
on fpi i the fraction of the 0.5-10.0 unabsorbed flux due to the power-law component, were derived from the best-fit model for 
the normalization N p i held fixed its 90% confidence limit. 



Table 6. Joint Spectral Parameters (Observations 1-4) 















fpl 




Obs. # 


Nh,22 


(eV) 


(km/(D/ 5 kpc)) 


a 


(phot/cm- 2 s _1 1 keV) 


(96) 


X? (prob) 



Observations 1-4 (0.6-9.0 keV); fcT c ff j00 varies 

1 0.42±0.03 121+.* 3 18A -3.t L1 ±0.6 (l-7tj;!)xl0- 5 13 1.0/41 (0.40) 

2 " 105tg " " " 23 
3+4 " 110+" " " " 19 

Observations 1-4 (0.6-9.0 keV); Ah varies 

1 n/a n/a n/a n/a n/a n/a 3.0/41 (2xlO -10 ) 

2 n/a 
3+4 n/a 

Observations 1-4 (0.6-9.0 keV); R OO Vetri6S 

1 0.39+J 122+j 3 17 - e -tl °- 7 ~o.t (0.8tJ;4)xlO" 6 13 1.13/41 (0.26) 

2 " " 12 -5l?'8 " " 23 
3+4 " " 14 -3±2l " " 18 

Observations 1-4 (0.6-9.0 keV); i*x,p!varies 

1 n/a n/a n/a n/a n/a n/a 1.60/41 (0.0085) 

2 " " ' " n/a n/a 
3+4 " " " n/a n/a 

Observations 1-4 (0.6-9.0 keV); a and Fx pi vary 

1 0.52+0.04 88±]P 32.5^ 3.2+0.2 (4.5+0.5) x 10" 4 54 1.08/39 (0.34) 

2 " " " 1.0+0.7 (9tg 7 )xl0" 6 10 
3+4 " " " 2.0+0.3 (9.5+3) xl0~ 5 30 

Observations 1-4 (0.6-9.0 keV); Ah and a fixed, all other values vary 

1 (0.34) 142+. 4 H-9±a3 <1.6xl0" 5 <16 0.73/41 (0.89) 

2 " 122+ 4 " " <l.lxl0" 5 <20 
3+4 " 128+4 " " (2.0ig 4 )xl0^ 5 28^ 

Observations 1-4 (0.6-9.0 keV); kT e ff varies, Ah and Roo fixed 

1 (0.34) 138+1 (13.0) 0.6+0.4 (0) (0) 1.30/43 (0.11) 

2 " 120+1 " " (0) (0) 
3+4 " 124+| " " (l.l+J^)xl0 _s 29 

Observations 1-4 (0.3-9.0 keV); fcT e fj varies, A^h an d Roc fixed 

1 (0.34) n/a (13.0) n/a (0) (0) 2.63/51 (2xl0" 9 ) 

2 " n/a " " (0) (0) 
3+4 " n/a " " n/a n/a 

Observations 1-4 (0.3-9.0 keV); kT e g varies, TVh and i?oo float 

1 n/a n/a n/a n/a (0) (0) 1.62/49 (0.0038) 

2 " n/a " " (0) (0) 
3+4 " n/a " " n/a n/a 

a Observations 1-4 (0.3-0.45, 0.6-9.0 keV); kT cS varies, and R x float 

1 0.42lJJ;°3 130 + g 15 - 9 t2:9 0.7+0.2 (0) (0) 1.26/48 (0.105) 

2 » ' 113±| " ' " (0) (0) 
3+4 " 118+^ " " (l-4+o.g)xlO- 5 26 

Observations 1-4 (0.3-9.0 keV); fcT e g varies, A^h and i?oo float 

1 n/a n/a n/a n/a (0) (0) 1.66/47 (0.0030) 

2 n/a n/a " " (0) (0) 
3+4 n/a n/a " n/a n/a 



Note. — a See Table ^ for luminosities derived from these values. b . The 0.45-0.6 keV band is exc luded from the data fi t 
with these models. Uncertainties are la, upper-limits are 90% confidence. Assumed distance is 5 kpc(Rutledge et al. 2001a). 
Notation "n/a" (not applicable) is used where spectral parameters were used in a fit which is not statistically acceptable. 
Upper-limits on f p i, the fraction of the 0.5-10.0 unabsorbed flux due to the power-law component, were derived from the 
best-fit model for the normalization N p i held fixed its 90% confidence limit. 



Table 7. X-ray Luminosities 



Obs. 


Thermal L x 


Total L x 


1 


4.6 


4.6 


2 


2.5 


2.5 


3+4 


3.0 


4.0 



Note. Assumed source 

distance d=5 kpc. Luminosities 
are corrected for absorption, in 
units of 10 33 ergs" 1 , 0.5-10 keV. 



Table 8. Observed vs. Best-Fit Model Counts in 0.45-0.6 keV 

Bin 



Obs.(i) 


Observed Counts (xj) 


Model Counts (/ij) 


1 


23 


26.0 


2 


23 


19.5 


3 


13 


21.4 


4 


10 


26.7 



Note. — Best-Fit Model parameters used to derive these values 
are in Table ^ the best-fit model for observations 1-4, using 0.3- 
9.0 keV, ignoring 0.45-0.6 keV, permitting A;T e fj j00 to vary between 
observations, and A% and -Roo to be the same between observations. 
Best fit ^=1-26/48 dof, prob=0.105. 



